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ABSTRACT 



Context. In a galaxy cluster, the evolution of spiral galaxies depends on their cluster environment. Ram pressure due to the rapid 
motion of a spiral galaxy within the hot intracluster medium removes the galaxy's interstellar medium from the outer disk. Once the 
gas has left the disk, star formation stops. The passive evolution of the stellar populations should be detectable in optical spectroscopy 
and multi-wavelength photometry. 

Aims. The goal of our study is to recover the stripping age of the Virgo spiral galaxy NGC 4388, i.e. the time elapsed since the halt 
of star formation in the outer galactic disk using a combined analysis of optical spectra and photometry. 

Methods. We performed VLT FORS2 long-slit spectroscopy of the inner star-forming and outer gas-free disk of NGC 4388. We 
developed a non-parametric inversion tool that allows us to reconstruct the star formation history of a galaxy from spectroscopy 
and photometry. The tool was tested on a series of mock data using Monte Carlo simulations. The results from the non-parametric 
inversion were refined by applying a parametric inversion method. 

Results. The star formation history of the unperturbed galactic disk is flat. The non-parametric method yields a rapid decline of 
star formation i 200 Myr ago in the outer disk. The parametric method is not able to distinguish between an instantaneous and a 
long-lasting star formation truncation. The time since the star formation has dropped by a factor of two from its pre-stripping value is 
190 ± 30 Myr. 

Conclusions. We are able to give a precise stripping age that is consistent with revised dynamical models. 
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1. Introduction 

Depending on the environment in which they move, spiral galax- 
ies undergo different processe s that can mod ify their struc- 
ture significantly (see Boselli & Gavazzil l2006l and references 
therein): 

- gravitational eff^ects (e.g. tidal interactions in galaxy-galaxy 
encounters), 

- hydrodynamical effects (e.g. ram pressure stripping or ther- 
mal evaporation), 

- hybrid processes, i.e. those involving both types of eff'ects, 
such as preprocessing and starvation. 

The closest cluster of galaxie s in the northern he misphere is the 
Virgo cluster (d ^ 16.7 Mpc. lYasuda et al.lll997l) with a mass of 
M - 1.2 X IO'^Mq and a radius of about 2.2 Mpc dpouque et al] 
l200lh . Vii- go is an evolving cluster that is still dynamically ac- 
tive. The cluster-core is centered on M87, the most massive ellip- 
tical galaxy. Other subclumps are falling into the potential well 
of the cluster 

One im portant characteristic of Virgo s piral galaxies is thei r 
lack of gas jGiovanelli & Havnes 1983, C hamaraux et al.lll980h . 
The amount of atomic gas in Virgo spirals is up to 80% less 



than for field galaxies of the same size an d morphological type. 
Virgo sp irals show trun cated HI disks ( Giovanelli & HavnesI 
1983, Ca vatteet al.l 19901) with respect to their optical disks. The 
galaxies on radial orbits are on average m ore HI deficient than 
the ones on circular orbits ( Dressier .1986) . IChung et all (l2009h 
find long HI tails associated with spiral galaxies located at dis- 
tances from 0.5 to 1 Mpc from the cluster center For these cases, 
ram pressure stripping is the most probable cause. The interstel- 
lar medium (ISM) of a spiral galaxy that is moving inside the 
potential well of a cluster, undergoes pressure due to the intra- 
cluster medium (ICM) that is hot (ricM ~ 10^ - 10** K) and ten- 
uous ( picM ~ 10"'' - 10""* atoms cm"-'). If this pressure is higher 
than the restoring force due to the galactic potential, the galaxy 
los es gas from the outer disk. Quantitatively this is expressed by 
the lGunn&Gotd([T97l criterion: 



PlCMV„ai ^ 27rG2star2gas, 



(1) 



where picM is the density of the ICM, Vgai the peculiar ve- 
locity of the galaxy inside the cluster, and Ejtar and Egas are 
the surface density of stars and gas, respectively. Ra m pres- 
sure s tripping has been studied theore tically (e.g. VoUmer et al 
ISchulz & Struck 2001, Ouili s et al.i l2000riAbadi et al 
iRoediger & Briiggeni j2008.) and observationally (e.g. 



2001 



1999 
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Kennev et alJ l2004 ISolanes et alJ l200lL ICavatte et alJ Il990i 
Warmelslll988ilChung et al.ll2009h . 
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NGC 4388 is a highly incUned Seyfert 2 spiral galaxy (of 
type Sab) with niB = 12.2 mag and a radial velocity of V^d ~ 
1400 km s"' with respect to the cluster mean. It is located at a 
projected distance of 1 .3° (~ 400 kpc at 16.7 Mpc distance) from 
the Virgo cluster ce nter (M87). NGC 4388 has lost 85% of its HI 
gas mass (ICavatte et al. 1990). The HI distribution is strongly 
tru ncated within the optic al disk. By observing the galaxy in 
Ha IVeilleux"etal] (Il999bh found a large plume of ionized gas 
extending 4 kpc above the pla ne. In subsequent SUBARU ob- 
servations lYoshida et all (12002) reveal a more extended tail out 
to 35 kpc to the northeast. iVoUmer & Huchtmeier (2003) per- 
formed 21 -cm line observations with the Eifelsberg 100-m ra- 
dio telescope. They discovered neutral gas associated with the 
Ha plume out to at least 20 kpc NE of NGC 4388's disk, with 
an HI mass of 6 x 10^ Mp. W ith interferometic observations 
lOosterloo & van GorkomI (l2005h show that this HI tail is even 
more extended, with a size of 1 10x25 kpc and a mass of 3.4x 10*^ 
Mo. 

In the case of NGC 4388. IVollmer & Huchtm"eie3 ( |2003b es- 
timate that ram pressure stripping is able to remove more than 
80% of the g a laxy's ISM, consistently with the observations of 
ICavatte et alj ( Il990l) . They also estimate that the galaxy passed 
the cluster center ~ 120 M yr ago. However this e stimation is 
based on the observations of lVoUmer & Huchtmeieri(,2003). The 
more extended HI tail found bv lOosterloo & van GorkomI il200l 
implies that this timescale increases. 

Once the gas has left the galactic disk, star formation, which 
is ultimately fueled by the neutral hydrogen, stops. Stripping 
is believed to progress inwards from the outermost disk on 
timescales of 10^ years, but at any given radius, stripping hap- 
pens on a shorter timescale. The stellar populations then evolve 
passively, and this is detectable both in optical spectra and in the 
photometry. The detailed analysis of the stellar light provides 
essential information on the stripping age, i.e. the time elapsed 
since the halt of star formation. These constraints are indepen- 
dent of dynamical models and can be compared to those de- 
rived from the gas morphology and kinematics. Crowl & Kenney 
(2008) analyze the stellar populations of the gas-free outer disk 
of NGC 4388, 5.5 kpc off the center (with our adopted distance 
of 16.7 Mpc). They used SparsePak spectroscopy and GALEX 
photometry. Their age diagnostics are based on Lick indices and 
UV fluxes, which they compare to stellar population synthesis 
predictions (models of Martins et al. 2005). Assuming that the 
past star formation rate in the galaxy has been constant on long 
timescales, they conclude the stripping of the gas from the outer 
disk occurs 225 + 100 Myr ago. Our renewed analysis is based on 
VLT FORS2 long-slit spectroscopy of NGC 4388 taken at two 
positions: 

- The first slit was pointed at 1.5 kpc from the center, in a gas- 
normal star-forming region of the disk. We did not take the 
exact center to avoid bulge contamination. 

- The second slit was pointed at 4.5 kpc from the center, just 
outside the star-forming gas disk. 

The spectra were combined with multi-wavelength photometry. 
The long-term star formation history of the galaxy was obtained 
using an extension of the non parametric inversion method of 
Ocvirk et al. (2006). It combined the information provided by 
the spectroscopic and photometric data to determine the star for- 
mation history using minimal constraints on the solutions. The 
results were then refined using a parametric method that assumes 




..J 



Fig. 1. B-band POSSII image of NGC4388 with HI contours. 
The sht positions for the inner and the outer regions are overlaid. 



the halt of star formation at the observed outer position in the 
disk occurs quasi-instantaneously. 

The paper is structured as follows. In Sect.s|2]and[3]we de- 
scribe the observations and explain how we extract the photom- 
etry that we use as input data. In Sect. |4] we introduce the new 
approach used in this paper that combines spectral and photo- 
metric analysis. In Sect. |5]and|6] we explain the results of this 
new method in the case of NGC 4388. Finally, in Sect.Qwe give 
our conclusions and compare our results to previous work. 



2. Observations 

2.1. Data set 

We performed observations of NGC4388 on May 2, 2006 at the 
European Southern Observatory Very Large Telescope (VLT) 
facility on Cerro Paranal, Chile. The instrument we used was 
the FOcal Reducer and l ow dispersion Spectrograph 2 (FORS2, 
Ap penzeller et al.l 11998 ) in long-slit spectroscopy observing 
mode. The detector system consisted of two 4096 x 2048 CCD. 
We chose 2x2 binning of the pixels (image-scale 0.252"/ 
binned pixel) and high-gain readout. We selected grism GRIS- 
600B with a wavelength range of 3350-6330 A and a resolution 
of 1.48 A / binned pixel. The data were acquired through a 1" 
slit, yielding a resolving power of A/AA » 780 at the central 
wavelength. 

The two slit positions used for the acquisition are shown in 
Fig.[T]together with HI contours. Individual exposure times were 
600 s for the inner region and 1350 s for the outer. The total in- 
tegration time for each region were 1800 s and 9450 s (Table 
[TJ. We used LTT 4816 and LTT 7987, two white dwarfs with 
my — 13.96 and my — 12.23, respectively dBakos et al.ll2002l 
iPerrvman et al.l[T997l) as spectrophotometric standards. 
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Table 1. Journal of observations (ESO program ID 77.B-0039(A)). Seeing, airmass, and UT (Universal Time) refer to the beginning 
of each acquisition. 



Observation 


Frame 


Exp. Time (sec) 


Seeing 


Airmass 


Date 




UT 


LTT4816Ph.St. 


1 


35 


0.91" 


1.484 


2006 May 


2 


23: 


:06 


NGC4388 Inner 


1 


600 


0.9" 


1.688 


2006 May 


2 


23: 


:45 


NGC4388 Inner 


2 


600 


1.11" 


1.616 


2006 May 


2 


23: 


:56 


NGC4388 Inner 


3 


600 


1.10" 


1.554 


2006 May 


3 


00: 


:07 


NGC4388 Outer 


1 


1350 


0.89" 


1.330 


2006 May 


3 


01 


:11 


NGC4388 Outer 


2 


1350 


0.83" 


1.289 


2006 May 


3 


01 


:34 


NGC4388 Outer 


3 


1350 


0.96" 


1.264 


2006 May 


3 


01 


:58 


NGC4388 Outer 


4 


1350 


0.67" 


1.256 


2006 May 


3 


02: 


:21 


NGC4388 Outer 


5 


1350 


0.61" 


1.423 


2006 May 


3 


04: 


:10 


NGC4388 Outer 


6 


1350 


0.55" 


1.519 


2006 May 


3 


04: 


:33 


NGC4388 Outer 


7 


1350 


0.55" 


1.651 


2006 May 


3 


04: 


:57 


LTT7987 Ph.St. 


1 


40 


1.34" 


1.010 


2006 May 


3 


10: 


:35 



2.2. Data reduction 

Data were reduced with the images and noao . twodspec pack- 
ages of the Image Reduction and Analysis Facility (IRAF13 
We created average bias and flat-field images us- 
ing daytime calibration, and we removed cosmic ray hits using a 
3cr-clipping rejection. 

After bias and flat-field correction, we combined the individ- 
ual galaxy exposures and applied the wavelength calibration to 
the spectral axis of the images. The aperture used for extracting 
the ID spectra covers a length of 32.76" along the slit for both 
the inner and the outer regions of NGC 4388. The sky subtrac- 
tion was based on a linear fit to the sky windows on either side 
of the aperture. Fin ally, we used t he spe ctrum of LTT 4816 and 
its intrinsic fluxes (iHamuv et al.L 1 19921) to calibrate the flux of 
the spectra. The flux calibration is ^ 20% accurate, based on 
the comparison of the energy distributions obtained using either 
LTT4816 or LTT7987 as a standard. 

The final reduced spectra for the inner and outer regions are 
shown in Sect. |5] The average signal-to-noise ratio per pixel for 
the inner and outer regions are x 60 and ^ 26, respectively, 
inside the wavelength ranges used in the analysis. 



3. Photometry 

For the photometry we used the following archive data: 



- GALEX (Gil de Paz et al. "2007"): FUV and NUV. 

- SDSS (Rel. 6, Adelinaii-McCarthy et al. 2008): u',g',r',i', 
andz'. 

- 2MASS dSkrutskie et al.ll2006h : J, H, and K. 

The archive images were resampled onto the coordinate grid 
used in F0RS2 observations using ALADIN (Bonnarel et al. 
[2000 ). The photometric zero points were derived using the star 
2MASS 7122549.86 + 12404 7.9 (Table|2]l. Fo r the G ALEX im- 
ages, we used the formulae of iMorrissev et al.l (l2007h to convert 
counts per second (CPS ) into AB magnitu des. The defined zero 
points are accurate to within +10% (Morrissey et al.ll20()7|). 

The exact locations of the two spectroscopic apertures on the 
images were determined by comparing the wavelength-averaged 
profile along the FORS2 slits with cuts through the SDSS g' im- 
age. In both cases a clean peak in the cross correlation function 



' IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 



for Research in Astronomy, Inc.. 
National Science Foundation. 



under cooperative agreement with the 



allowed us to identify the slit position to within 1 pixel size. The 
GALEX images have a spatial resolution of about 4". Light mea- 
sured within a 1 " aperture is contaminated by neighboring areas, 
but this effect is less than the unceitainties already accounted for, 
which are large because of low signal levels. 

We estimated the uncertainties on the final magnitudes for 
the galaxy apertures from a combination of the uncertainties 
on the flux measurements of the reference star (zero point) and 
those of the positioning of the slit (see Table|2l): 



cr. 



TOT 



SLIT 



+ cr; 



STAR- 



(2) 



4. Method 

The spectral energy distribution (SED) of a galaxy can be con- 
sidered as the integrated light produced by different stellar pop- 
ulations with different ages and metallicities. The light and mass 
contributions of each population depend on the star formation 
history of the observed region during the galaxy's life. 

If we define an initial mass function and a set of isochrones, 
we can obtain the intrinsic spectrum B°(/l, t, Z) of the single stel- 
lar population of age t, metallicity Z, and unit mass by integrat- 
ing over the stellar masses. Assuming that the metallicities of the 
stars can be described by a single-valued age-metallicity relation 
Z(f), we can derive the unobscured SED of a galaxy at rest: 



^'rest(^) 



*^'niin 



SFRit) B\A,t,Z{t)) dt. 



(3) 



in which S FR(t) represents the mass of new stars born per unit 
time, with the convention that fmi,, = is today and fmax the 
Hubble time. 

Since we observe the light and not the mass in a galaxy, it 
is more convenient to convert the mass-weighted spectral ba- 
sis B"(/1, f, Z) into a luminosity-weighted basis. The difference is 
that B"(/1, t, Z) defines the spectrum of a single stellar population 
(SSP) of unit mass, and B(A,t,Z) defines an SSP spectrum of 
unit flux. Instead of mass contributions we deal with light con- 
tributions, thus converting the SFR into the luminosity-weighted 
stellar age distribution (hereafter SAD): 



SFR(t) 
AA 



Z)dA, 



(4) 



4 
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Table 2. NGC 4388 photometry. 1 " - 2""' row: magnitude and eiTor for the reference star 2MASS 7122549.86 + 124047.9. 3''' - 4"' 
row: magnitude and error for the inner region. 5''' - 6''' row : magnitude and error for the outer regions. 7''' - 8''' row: total uncertainty 
for inner and outer regions. For the SDSS and 2MASS filters the magnitude are given in AB and Vega systems, respectively. 



Filter 


FUV 


NUV 


u' 


8' 


r' 


i' 


z' 


J 


H 


K 


m 






16.78 


15.06 


14.35 


14.69 


14.07 


13.08 


12.65 


12.52 








0.025 


0.034 


0.017 


0.002 


0.016 


0.021 


0.03 


0.024 


Inner 


20.40 


19.95 


18.45 


17.10 


16.37 


16.03 


15.76 


14.38 


13.60 


13.32 


(Tsur-Inner 


0.17 


0.04 


0.01 


0.02 


0.04 


0.04 


0.04 


0.04 


0.02 


0.05 


Outer 


22.21 


21.30 


19.49 


18.07 


17.47 


17.18 


17.03 


16.08 


15.06 


14.75 


o-iur-Outer 


0.31 


0.07 


0.08 


0.06 


0.07 


0.07 


0.03 


0.33 


0.04 


0.12 


cr 7-or -Inner 


0.2 


0.1 


0.03 


0.04 


0.04 


0.04 


0.04 


0.05 


0.04 


0.06 


(Tror -Outer 


0.33 


0.12 


0.07 


0.07 


0.07 


0.07 


0.03 


0.33 


0.05 


0.12 



where AA - /l,nax - ^min is the available wavelength domain, and 
A(f) gives the contribution to the total light from the stars of age 
[f,f-i-df]. With 



B(A,t,Z) 



B"(A, t,Z) 



ijj;-BO(A,t,Z)dA 



(5) 



Eq.[3]becomes 

FrM = J A{t)B(A, t, Z(t)) dt. 



(6) 



In the same way we can associate a photometric value with each 
spectrum B(A, t, Z(t)) for the set of bandpasses defined in Table 

m 



BphotCy, t,z{t)) = 



J^°'" B(A,t,Z(t)) ■ Ty(A) A dA 



T,(A) A dA 



(7) 



in which y = [FUV,NUV,u',g',i',r',z',J,H,K], and T,. is the 
transmission curve associated with each y. Unobscured photom- 
etry Fphotiy) is given by 



= j'"'"A(f)Bphot(3',f,Z(f)) df. 



(8) 



Eq.|6]is not complete, because it does not take other effects into 
account that can modify the final shape of the SED of a galaxy, 
as follows: 

• Extinction: when fitting spectroscopic data we use a flexi- 
ble continuum correction that can account both for the reddening 
due to dust and for flux calibration errors. The method allows us 
to use the information present in the spectral lines without us- 
ing the continuum of the spectrum, which is preferable when 
the flux calibration is not perfect. In practice, we define a set of 
equally spaced anchor points across the wavelength range of the 
optical spectra. Their ordinates are then optimized in such a way 
that the cubic spline interpolation through the points cancels out 
any SED difference between models and data. We refer to this 
adjustable correction as NPEC hereafter (for non parametric es- 
timate of the continuum). 

For the photometry, however, extinction mu st be accounted 
for explicitly. We adopted the attenuation law of lCalzettil (1200 ih 
that uses the colour excess E(B - V) as a single parameter (see, 
however, Sect.|6]l. 

• Radial velocities of the stars: we assumed that the veloci- 
ties of stars of all ages along the line of sight have the same ve- 
locity distribution. We can approximate the spectrum of a galaxy 



as the convolution of the spectrum at rest and a line of sight ve- 
locity distribution g(v): 



f"™' I A \ dv 

f rest , ^ ■ ]-8(v) , 

J,, , V 1 H- v/cJ 1 + v/c 



(9) 



where Frest is defined in Eq. |6] and c is the velocity of light. The 
integrals in Eqs.|6][8] and |9] are transformed into discrete sums 
for numerical treatment. 

Recovering the star formation history ultimately means find- 
ing a A(f) that fulfills Eqs. |6] and |8] while accounting for the 
corrections just described. To do this, we used two different ap- 
proaches: 

- a non parametric method in which we recover the star for- 
mation history by resolving t he associated invers e problem 
with regularized methods (see lOcvirk et al.ll2006al) : 

- a parametric method in which we define a set of possible A(f) 
depending on one free parameter, the time elapsed since the 
stripping event. We then find the most probable solution by 
minimizing the classical function. 

The non parametric method has the advantage of providing the 
star formation history and trends in the age-metallicity relation 
of the galaxy with minimal constraints on their shape. On the 
other hand, the regularization of the problem did not allow us to 
recover functional forms with large gradients, such as those ex- 
pected for a ram pressure stripping event. This is why we com- 
bined the results from the non-parametric analysis with a para- 
metric analysis. 

We chose as SSP library the models of iBruzual & Chariot! 
(I2003h that cover a time interval Af = [0 - 19.5] Gyr, 
a wavelength range AA - [100 - 24850]A with an aver- 
age spectral resolution (FWHM) R - 2000 at optical wave- 
lengt hs. These SSP spectra a re constructed with stellar spectra 
from iLe Borgne et aL (12003 ). The unde r lying stellar evolutio n 
tracks are those of Alongi et al. ('1993), Bressan et al. ( 1993|, 
Fagotto et al- (.1994a) . .Fagotto et al.. (.1994b) and .Girardi et all 
(I1996h . 



4. 1 . Non parametric method 

The non para metric method used is a modificatio n of 
STECKMAlfl (lOcvirk et all l2006aL lOcvirk et all l2006bl) . ex- 
tended to deal with spectroscopy and photometry jointly. 
Assuming gaussian noise in the data, we estimate the most likely 



http://astro.u-strasbg.fr/~ocvirk/ 
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solution by minimizing the following 2^(X) function for the star 
formation history and age-metallicity relation: 

Q,{X) = (1 - a) ■;rspec(X) + a ■ A„,(X) ■ TO, (10) 



in which 



X is a vector that includes the SAD, the metallicity evolu- 
tion, the color excess E{B— V) for the photometry, the NPEC 
correction vector for spectral analysis and a spectral broad- 
ening function (BF). The last includes the effects of the line 
of sight velocity distribution and the broadening introduced 
by the instrumentation. 

4ec(X)=^i:^ ^, (11) 

is the associated with the spectrum. Fmdispec('i) is the 
model with a defined X and F^^^^iX) is the observed spec- 
trum, <Tspj.(, is the associated error in the observations, and 



A'^ ~ 1700 is the number of points in the spectrum, which 
we can consider approximately equal to the number of de- 
grees of freedom of the problem. 

1 (■Fmdlphot(3') - ■P'photC}')) 

;^phot(X) = ^ 2 ^ -2 (12) 

y i=\ "phot 

is the;if^ associated with the photometry, where FmdiphotCy) is 
the model photometry obtained from Eq.[8]with a chosen X, 
Fptiotiy) are the photometric data and cr^^^^ are the errors in 
the measures, and A^y represents the number of bandpasses 
used. 

^lP{X) = ^JL,P{SAD) + jUzf'(AMR) (13) 

is a pena lty function needed t o regularize the problem. As 
shown in Ocvirk et al. (2006a), the inverse problem associ- 
ated with Eqs.|6]and[8]is ill-posed. The function P(X) is cho- 
sen to yield high values when the SAD and AMR are a very 
irregular function of time or when the BF or the NPEC cor- 
rection are too chaotic. The set of /i = (/ijc/iz./^wPEC./^BF) 
are adjustable parameters that control the weight of each 
f (X) in the final estimation. 

a determines the relative weights of the photometric and 
spectroscopic constraints. 



4.2. Parametric method 

In this case we assume an exponential star formation history 
before the stripping event, a metallicity and a BF, all consis- 
tent with the non parametric results. We reproduce the strip- 
ping by a linear decrease in the SFR on timescales between 
< f < 100 Myr. For each of these truncated star formation 
histories, we calculate the NPEC correction and E(B - V) that 
produce the lowest The minimum of x^{t) is taken to pro- 
vide the most likely stripping age, which we define as the time 
elapsed since star formation has dropped by a factor of two from 
its pre-stripping value. The spectroscopic and photometric con- 
tributions to the are weighted as before: 



xiiiX, f) = (1 - a) ■;rLc(X, t) + a-x^ 



phot 



(X,f), 



(14) 



where a,xlpec(^' 0; andXpiiot(-X' have the same meaning as in 
Sect.OI 



5. Results 

5.1. Non parametric Inversion 

First we point out that the weightings fj. for the penalty are a cen- 
tral issue of the non parametric method and that their determi- 
nation is not a trivial problem. T here are different ways of fixing 
these values dTitterin gtonLll985h . In the case o f nonlinear prob- 
lems, one has to proceed via empirical tuning dCraig & Brownl 
1986) to define the values of fi below which the results present 
artifacts, and above which the smoothness of the solution is com- 
pletely due to the penalty. 

We chose as penalty f unctions two finite d iff'erence opera- 
tors (L - Di and h - Dz inl Qcvirk et al.ll20 06a) that favor small 
first derivatives of the AMR and small second derivatives of the 
SAD and BF. The penalties for the star formation history and the 
metallicity evolution are defined in logarithmic timescale. We 
explored the effect of the weight coefficients /i through a cam- 
paign of inversions of artificial spectra. For the NPEC correction, 
we found it was not necessary to require smoothness explicitely, 
and the penalization simply acts to normalize the continuum cor- 
rection (thus avoiding the degeneracy between the absolute val- 
ues of this correction and the star formation rates). The method 
recovers similar well-behaved BF for a wide range of fisF- The 
campaign showed that it is not possible to recover more than a 
tentative linear trend in the AMR with the data available to us, 
and we chose to simplify the problem with constant metallicity 
(i.e. large fxz)- A realistic metallicity evolution shows a rapid in- 
crease at long lookback times an d little increase over the past 
5 Gyr (iBoissier & Prantzosll2000h . From our method we derived 
a time-bin ave raged metallicity over a Hu bble time (Z). For typ- 
ical models of lBoissier & PrantzosI (12000'). this mean metallicity 
is AZ ~ 0.005-0.01 lower than the time-bin averaged metallicity 
over the past 5 Gyr 

For the stellar age distribution, we adopt the smallest jj ^ that 
provides robust solutions, i.e. low sensitivity to the noise in the 
data. This choice was based on successive inversions of pseudo- 
data with different /i j . The stability of the results has been tested 
via Monte Carlo simulations: we added a Gaussian noise with cr 
equal to the noise of the data to the input spectrum and repeated 
the inversion. In this way we fixed the limits for the fi^ value. 
In the campaign of tests, we also verified that, within reasonable 
limits, the shape of the initial guess does not affect the recov- 
ered solut ion significantly. We ended up taking a semi-analytical 
model of Boissier & PrantzosI (l2000h as an initial guess for the 
star formation rate and the age-metallicity relation, and a con- 
stant for the BF and the NPEC correction. 

For each region we reconstructed the star formation history 
using the VLT spectrum alone (a - 0), the photometry alone 
(ff = 1) and spectrum/photometry at the same time (a - 0.5). In 
the combined analysis the choice of a = 0.5 gives, in the final 
estimation of 2(X), the same weight to the spectral and photo- 
metric analysis. Our conclusions are not significantly affected by 
the choice of a, as verified through a campaign of inversions. 

5.1 .1 . Inner Region 

- VLT spectrum (a - Q). 

The spectrum shows a fit (top panel in Fig. O with xlpec - 
2.2 and a star formation history (Fig. |3]l that is approxi- 
mately flat. The time-bin averaged metallicity is nearly solar, 
(Z) = 0.018 + 0.003. For comparison we took three emis- 
sion fines ([OII]3727, Hy64861, [OIII1500 7) to recove r the 
metallicity using the strong line method (lPilvuginl2000l) . We 
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Fig. 2. NGC 4388 inner region fits. Top panel: spectra (black line) and best fit (red line) obtained from inversion with a - 0. The 
bottom of the panel shows the residuals (magenta line) and the cr (green line). The cyan vertical lines show the mask used for 
emission and sky lines. Middle panel: photometry (black crosses) and best fit (red crosses) obtained from inversion with c = 1 
overplotted to the corresponding flux (green dotted line). In the bottom of the panel are shown the residuals (magenta crosses). 
Bottom panel: same as the middle panel, but using a - 0.5. We do not show the fit of the spectrum for the case a = 0.5, because it 
is indistinguishable by eye from the top panel. 
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Fig. 3. Non parametric spectral inversion (a - 0) of the inner 
region of NGC 4388. Top and middle panel: star formation 
history (SFH) and luminosity weighted stellar age distribution 
(LWSAD) vs. lookback time. Bottom panel: spectral broadening 
function (BF). In each panel the solid lines shows the results of 
the inversion with associated error bars obtained from 30 Monte 
Carlo simulations. The dashed line shows a constant star forma- 
tion rate. All the figures are in arbitrary units. 
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Fig. 4. Non parametric photometric inversion (a = 1 ) of the in- 
ner region of NGC 4388. Top and bottom panel: star forma- 
tion history (SFH) and luminosity-weighted stellar age distribu- 
tion (LWSAD) vs. lookback time. For a detailed description see 
Fig.S 



obtained (Z) x 0.021 consistent with the results of the inver- 
sion. The broadening function is centered on -100 km/s. The 
shift and the width of the broadening function are consistent 
with expectation based on gal actic rotation and non -circular 
motions in a barred potential (IVeilleux et al.lll999bh . 
Photometry (a = 1 ) 

The fit reproduces the observations well with a = 0.36 
(middle panel of Fig.|2]). The star formation history (Fig.|4]i 



Fig. 5. Non parametric combined inversion (a - 0.5) of the in- 
ner region of NGC 4388. Top and middle panels: star formation 
history (SFH) and luminosity-weighted stellar age distribution 
(LWSAD) vs. lookback time. Bottom panel: spectral broadening 
function (BF). For a detailed description see Fig. [3] 



is quite similar to the star formation history recovered in the 
case a -Q, except in the past 10 Myr. The metallicity is con- 
sistent with the spectroscopic results with larger error bars. 
From the photometric inversion we recovered the average 
reddening of the stars in the region, E{B - V) = 0.18. 
VLT spectrum + photometry (a = 0.5) 
The fit of the spectrum is essentially identical to the case 
a = 0, but we do not show it. The photometry is shown in the 
bottom panel of Fig. |2l The total x^oi - 1-4 is intermediate 
between the two former cases. The star formation history and 
the luminosity-weighted stellar age distribution are shown in 
Fig. 13 



5.1.2. Outer region 

- VLT spectrum (a - 0) 

We obtain a fit (top panel of Fig. |6j with a x%ec - 0.43. 
The recovered star formation history (Fig. |7]i is flat until a 
lookback time of ~ 200 Myr at most, then it shows a drop, 
as expected from the ram pressure stripping scenario. The 
time-bin averaged metallicity is (Z) = 0.01 ± 0.006, about 
half of that of the inner region. The broadening function is 
centered on ^ 200 km /s, consistent with galactic rotation 
(IVeilleux et al.l[T9"99bh . Its width is dominated by the spec- 
tral broadening of the instrument. 

- Photometry (a = 1 ) 

The star formation history shows a shape similar to the a - 
case, with a departure from a constant value at a lookback 
time of ~ 200 Myr at most (Fig.[8]l. The extinction is E{B - 
V) = 0.07 and the ;fphot = 1-1 (middle panel of Fig.|6]l. 

- VLT spectrum + Photometry (a - 0.5) 

The total fit (bottom panel of Fig.|6]l has aXtot - 0-51. The 
spectral fit is very similar to the case = so is not shown. 
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Fig. 7. Non parametric spectral inversion (a - 0) of the outer 
region of NGC 4388. Top and middle panels: star formation 
history (SFH) and luminosity-weighted stellar age distribution 
(LWSAD) vs. lookback time. Bottom panel: spectral broadening 
function (BF). For a detailed description see Fig. [3] 

The total star formation history was flat until ~ 200 Myr 
ago, then it decreases steeply (Fig. |9]i. The gas stripping of 
NGC 4388 truncated the star formation history between 100 
and 500 Myr ago. 

Fig EHcompare the recovered star formation history for the 
inner and the outer regions to a flat star formation. While the 
inner region shows good agreement, the outer region, some- 
where between 100-200 Myr, significantly deviates from a 
flat star formation history. Since we regularize the inverse 
problem by requiring the solution to be smooth, instanta- 
neous truncation is rejected by construction. Therefore we 
cannot assess whether the progressive nature of the decline 
of the post-stripping star formation is real or if it is a con- 
sequence of the penalization. We get back to this point in 
Sect.lO 



As expected, the non parametric method 

1 . provides constraints on the long-term star formation history 
of the galaxy, 

2. confirms a recent radical change in the star formation history 
of the outer disk, 

3. does not provide a sharply truncated star formation history. 



5.2. Parametric inversion 

We applied the parametric method to the spectrum of the outer 
region of NGC 4388 to quantify the time elapsed since star for- 
mation dropped by a factor of two from its pre-stripping value, 
which we call stripping age. For the truncation of star forma- 
tion, we used star formation rates that decline linearly during 
to 100 Myr. For the pre-stripping star formation history and 
metallicity, we used the results of Sect. 15.11 This is essential be- 
cause, as explained below, the derived stripping age depends on 
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Fig. 8. Non parametric photometric inversion (a - 1) of the 
outer region of NGC 4388. Top and bottom panels: star forma- 
tion history (SFH) and luminosity-weighted stellar age distribu- 
tion (LWSAD) vs. lookback time. For a detailed description see 
Fig.E 
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Fig. 9. Non parametric combined inversion (a - 0.5) of the outer 
region of NGC 4388. Top and middle panels: star formation 
history (SFH) and luminosity-weighted stellar age distribution 
(LWSAD) vs. lookback time. Bottom panel: spectral broadening 
function (BF). For a detailed description see Fig. [3] 



the luminosity-weighted Z and on the ratio of young-to-old stars 
at the time of stripping. The first point to clarify is the impor- 
tance of the choice of metallicity evolution in the stripping age 
determination. 

To test the influence of different AMRs on determination of 
the stripping age, we used a grid of constant metallicity mod- 
els, as well as two AMRs from the galaxy evolution models of 
iBoissier & PrantzosI (l2000h . We find that the stripping age de- 
creases when the average metallicity increases. The stripping 
age of the parametric method is only sensitive to the metallic- 
ity averaged over the past 5 Gyr Since the time-bin averaged 
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Fig. 6. NGC 4388 outer region fits: top panel: spectra (black) and best fit (red) obtained from inversion with a = 0. The bottom 
of the panel shows the residuals (magenta) and the observational errors (green). The cyan vertical lines show the mask used for 
emission and sky lines. Middle panel: photometry of NGC 4388 (black crosses and error bars) and best fit (red crosses) obtained 
from inversion with a - \. The best model spectrum is overplotted. The photometric residuals are also shown (magenta crosses). 
Bottom panel: same as the middle panel, but using a - 0.5. We do not shown the fit of the spectrum for the case a - 0.5, because it 
is indistinguishable by eye from the top panel. 
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metallicity from the non-parametric method, which is averaged 
over a Hubble time, is lower by AZ ~ 0.005-0.01 (cf. Sect. ED, 
we took a constant metallicity of Z = 0.018 for the parametric 
method. We studied the;^'^ as a function of the stripping age. As 
done for the non parametric method, we applied the paramet- 
ric method to the VLT spectrum, the photometry alone and then 
spectrum and photometry together We find that the stripping age 
is independent of the duration of the star formation truncation. 
The analysis does not allow different durations of the star for- 
mation truncation to be distinguished. In the following we only 
discuss the case of instantaneous truncation. 

- VLT spectrum {a - 0) 

The top panel of Fig. [TO] shows the value qI}^ as a function 
of the stripping age, yielding a most likely stripping age of 
190 Myr 

- Photometry (a -Y) 

We obtained '^x''^^^.^ - 0.54 (middle panel of Fig. [TOt and a 
stripping age of 190 Myr consistent with the result obtained 
for the VLT spectrum analysis. We also recovered an E{B - 
y) = 0.1, in good agreement with the reddening obtained in 
the equivalent non parametric problem. 

- VLT spectrum + photometry {a - 0.5) 

The minimum is 0.57, and the stripping age is again of 
190 Myr (bottom panel of Fig. [TOt . 

With the non parametric method, we fixed an upper limit for 
the stripping age of about 100-200 Myr The parametric method 
improved the precision by giving a stripping age of 190 Myr It is 
remarkable that the photometric and spectroscopic data provide 
good agreement on the stripping age. Uncertainties on this value 
are discussed in the next section. 

6. Discussion 

6.1. Uncertainties in tine parametric mettiod 

In determining of the stripping age, a crucial point is to deter- 
mine the uncertainty of the result. To clarify this point, we dis- 
cuss here the influence of the potential sources of errors in the 
method. 

- Monte Carlo simulations 

We add a Gaussian noise to the spectral and photometric 
data, and perform 500 Monte Carlo simulations that way. 
As expected from statistics, the absolute values of the mini- 
mum vary and the dispersion is ~ 1 / Va^, where is the 
number of degrees of freedom (> 1700 for the spectroscopic 
analysis, 7 for the photometric analysis). However, the loca- 
tion of the minimum in the x^ curves varies very little. The 
error on the stripping age associated with the noise is 20 Myr. 

- Extinction law 

The extinction law of ICalzettil (1200 Ih was designed to 
provide a reasonable correction for starburst galaxies. 
Compared to other laws in the literature, such as the standard 
Milky Way extinction laws, it is relatively flat. The choice of 
the extinction law affects the UV-optical colors particularly 
strongly. However, in the external regions of NGC 4388 from 
which gas has mostly been removed by ram pressure strip- 
ping, the effect of the choice of the extinction law can only 
be small. As a test, we apply to the model spectra two differ- 
ent extinction laws (Calzetti 2001, CardelH et al. 1989), and 
we compare the recovered stripping ages. The results differ 
by 10 Myr. 
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Fig. 10. Parametric inversion of the outer region for an instanta- 
neous star formation truncation. Each panel shows the^^ evolu- 
tion as a function of the time elapsed since the truncation of star 
formation. To build the set of spectra, we used a constant star 
formation rate and metallicity (Z = 0.018). Top, middle, and 
bottom panels show the x^ evolution for spectral, photometric, 
and combined analysis, respectively. The bottom panel overplots 
the results of the spectral (dotted line) and photometric (dashed 
line) analysis. Labels indicate the underlying star formation his- 
tory, the metallicity, the values of the minimal and the result- 
ing stripping ages. 



- Long-term star formation history 

At fixed metallicity, we build families of model spectra with 
different star formation histories. We consider exponentially 
decreasing star formation rates (SFR oc e^'"/') with timescales 
T =6, 10, and 50 Gyr, as well as a flat SFR. The stripping 
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age increases from lOOMyr for r = 6Gyr to 190Myr for a 
flat SFR, because the ratio of young-to-old stars before the 
ram pressure stripping event is higher with flatter star for- 
mation rates. Shorter timescales are excluded, because the 
corresponding spectra are too heavily dominated by old stars 
even before any cut in the SFR is considered. In our study 
the long-term star formation history is constrained strongly 
by the non parametric inversion of Sect. l5.ll For a wide range 
of penalization weights, constant star formation rates are fa- 
vored. The uncertainty on the stripping age obtained when 
considering only regular star formation histories consistent 
with the non parametric results is reduced to about 10 Myr 
- Metallicity 

At fixed star formation, we build sets of model spectra using 
different metallicity values. The stripping age increases with 
increasing metallicity. With our study the time-bin averaged 
metallicity over the last 5 Gyr is established through the non- 
parametric analysis. The error on the metallicity is 20%. The 
resulting error on the stripping age is 10 Myr. It is worth not- 
ing that the agreement between the stripping ages determined 
from the photometry alone and from the spectroscopy alone 
is best with the flat star formation history and the quasi-solar 
metallicity found by the non-parametric analysis. With other 
assumptions, this agreement is lost. 

6.2. Gas dynamics and star formation 

The non parametric method fixed an upper limit for the strip- 
ping age of about 100-200 Myr, with a rather unconstrained du- 
ration for the truncation. Numerical simulations (Abadi et al. 
1999, Vollmer et al. 2001, Rodiger & Bruggen 2006) show that 
the ISM located at a given distance from the galaxy center is 
stripped rapidly within a few 10 Myr. Therefore in the fam- 
ily of possible interpretations of the outer disk data, we favor 
a posteriori the scenarios with sharp truncations, although the 
timescale of the decline cannot be constrained with the data pre- 
sented here. The instantaneous truncation scenario can only be 
studied with a parametric method, since the penalization used 
in the non-parametric method rejects luminosity weighted stel- 
lar age distributions with sharp variations. A rapid halt in star 
formation suggested by the dynamical models can only be stud- 
ied with a parametric method, where we make the simplifying 
assumption that the gas is stripped instantaneously. The derived 
stripping age of 190 Myr is consistent with the extent of the ob- 
served Hi tail (Oosterloo & van Gorkom 2005): a tail extent of 
80 kpc in the plane of the sky and along the line of sight with 
the given stripping age leads to a total velocity of ~ 570 kms"' 
and a radial velocity of ~ 400 kms"'. The latter corresponds 
to the difference between the galaxy's systemic velocity and the 
velocity of the Hi tail (Oosterloo & van Gorkom 2005). The full 
extent of the Hi tail was not known at the time when Vollmer 
& Huchtmeier (2003) presented their dynamical model. Their 
tail has an extent of 40 kpc with an associated stripping age of 
120 Myr. A new, revised dynamical model (Vollmer, in prep.) 
with the observed extent yields a ram pressure peak ~ 170 Myr 
ago, consistent with our findings. 

7. Conclusions 

VLT FORS2 spectroscopic observations of the inner star- 
forming and outer gas-free disk of the Virgo spiral galaxy NGC 
4388 have been presented. Previous observations indicate that 
this galaxy has undergone a recent ram pressure stripping event. 
Once the galaxy's ISM is stripped by ram pressure star formation 



stops. We detect this star formation truncation in the spectrum 
and the multiwavelength photometry of the outer disk region of 
NGC 4388. 

To derive star formation histories we extend ed the non para- 
metric inversion method of iOcvirk et al ] (l2006a ^ making a joint 
analysis of spectroscopic and photometric data possible. The 
new code was tested on a series of mock data using Monte Carlo 
simulations. We find that the results are stable, once that mini- 
mization has converged. The uncertainties for the young stellar 
ages (< 10 Myr) can be large because of the uncertainties of 
stellar models. Whitin reasonable limits, e.g. acceptable the 
shape and the normalization of the initial guess does not signif- 
icantly affect the recovered solution (at a fixed signal-to-noise 
ratio). 

The new inversion tool was applied to our spectroscopic and 
photometric data for NGC 4388. We explored the effect of differ- 
ent penalizations in case of spectral analysis, photometric analy- 
sis and combined analysis. The main results are (i) the recovered 
star formation history is flat for the inner disk spectrum and (ii) 
for the outer disk region the inversion yields a star formation 
drop at a lookback time of ~ 200 Myr at most. 

We have introduced a parametric method that refines the de- 
termination of the stripping age, the time elapsed since the star 
formation has dropped by a factor of two from its pre-stripping 
value. Based on the non parametric results, we assumed a flat star 
formation before the stripping event and an almost solar metal- 
licity. We approximated the effect of gas stripping by cutting 
or decreasing linearly the star formation at a different lookback 
time < f < 1 Gyr The obtained set of spectra was compared 
with the observed spectrum of the outer disk of NGC 4388. The 
effect of the potential sources of error in the stripping-age deter- 
mination were evaluated. 

The parametric method leads to a stripping age for NGC 
4388 of ~ 190 + 30 Myr. It cannot distinguish between an in- 
stantaneous and a longer lasting star formation truncation. The 
de rived stripping age agree s with the results of a previous work 
of ICrowl & KennevI (l2008l) and with revised dynamical models. 
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